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FIG. 1: Closeup of experimental fracture fronts from 

Alava and Zapperi (AZ) Q question whether the fracture fronts we observe in j2j are 
self affine. We use Family- Vicsek scaling to determine the two scaling exponents a and 
z. AZ claim that this is not enough to determine whether the front is self affine and they 
go on to point to the presence of overhangs as their evidence of fractality rather than self 
affinity In Fig. ^ we show details of experimental fracture fronts from Fig. 1 of jjj . There 
are significant overhangs, but these experimental fronts are self affine. In fact, as long as 
one allows a damage cloud to develop, overhangs are unavoidable. Family- Vicsek scaling 
implies non-isotropic scaling. A consequence is that the front width w scales with the width 
of the system L x . Non-isotropic scaling is the defining property of self-affine surfaces. 

AZ claim that the source of the /^.-dependency of the damage length scale l y , Eq. (9) in 
Q], comes from the rescaling of the model's elastic constants when changing L x . However, 
this rescaling is necessary to ensure that the elastic properties remain unchanged when the 
system size is changed under uniform loading conditions. It is caused by the non-local 
nature of the problem introduced by the Green function Gij, Eq. (4) in 2]. Under uniform 
loading condition, Eq. (5) in P reads Ui = u = J2j Gi,jfj — J2j Gijb 2 a. Here b is the lattice 
constant. If u, local deformation, and a, local stress, are to be independent of size, L x and 
L y , we must have that 

G it j = constant . (1) 
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Since Gij = Gi-j, we may for estetic reasons write J2j Gij — constant = Gij/ (L x x L y ), 
as there is no dependency on the index % in Eq. ((T} . This was the way we chose to present it 
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FIG. 2: Average stress (a) as a function of imposed displacement D from virgin system to complete 
failure. The rescaling of Gij, Eq. and of the threshold distribution ensures that the curves 
collapse for different system lengths L x and fixed L y = 128. 

in Q| . We point out again that both indices in the double sum run over L x x L y sites. 
We demonstrate the correctness of the rescaling in Fig. El where we show the collapse of the 
loading curves obtained for different system sizes after using Eq. (JTJ). Only when the elastic 
properties have been rescaled as just described, one may proceed to use finite size scaling as 
done in 0. 

We end this Reply by pointing out that Ramanathan and Fisher jj] measured numerically 
v = 1.52 ± 0.02 using a very different model, which by construction cannot produce fractal 
fracture lines. This is in complete agreement with our model, which gave v = 1.54. The 
fracture roughness exponents, however, are very different. 
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